Novel mechanism of photoinduced reversible phase transitions in molecule-based 

magnets 
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A novel microscopic mechanism of bi-directional structural changes is proposed for the photo- 
induced magnetic phase transition in Co-Fe Prussian blue analogues on the basis of ab initio quantum 
chemical cluster calculations. It is shown that the local potential energies of various spin states of 
Co are sensitive to the number of nearest neighbor Fe vacancies. As a result, the forward and 
backward structural changes are most readily initiated by excitation of different local regions by 
different photons. This mechanism suggests an effective strategy to realize photoinduced reversible 
phase transitions in a general system consisting of two local components. 
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Repeatable switching-on and -off of magnetization 
by external stimuli such as light is one of the most 
fascinating phenomena with potential applications in 
next generation's information storage and processing. 
A bi-directional photo-induced magnetization was first 
discovered in a cobalt-iron Prussian blue analogue, 
K .4Coi. 3 Fe(CN) 6 - 5H 2 0-||. Illumination of visible 
light (500 - 700 nm) at low temperature induces a bulk 
magnetization (presumably, ferrimagnetism) , which can 
be eliminated by illumination of near-IR light (~1300 
nm). In spite of various experimental and theoretical 
efforts (l]J^|| the microscopic mechanism of the re- 
versible magnetization is still not clear. In this Letter we 
report ab initio quantum chemical cluster calculations for 
the Co-Fe Prussian blue analogue, unveiling the micro- 
scopic mechanism of the bi-directional photo-induced lo- 
cal structural changes that trigger the phase transitions. 

The crys- 
tal of a Prussian blue analogue Ki_2 a; Coi+ a ;Fe(CN)6 is 
composed of two metallic sites located on vertices of the 
cubic lattice and each surrounded by six cyano moieties, 
as shown in Fig. [I]. The d-orbitals of transition metals 
split into ^2g and e g orbitals by the ligand field. In the 
case of x^0, there are vacancy sites with replacement 
of CN by H2O as shown in Fig. 0(b). Various fascinat- 
ing phenomena including a room temperature magnet 
H, electrochemically tunable magnets Jib]] , transparent 
and colored magnetic thin films [JQJ, and photo-induced 
magnetic dipole inversion [ p"2| have been observed in such 
non-stoichiometric compounds. We will show that this 
non-stoichiometric aspect is essential for the reversible 
photo-induced magnetization. 

The low spin (LS) configuration of the ground nonmag- 
netic state and the high spin (HS) configuration of the 
meta-stable magnetic state of K .4Coi.3Fe(CN) 6 - 5H 2 
are most likely Co 111 ^ 6 ,,? = 0)Fe n (de 6 , S = 0) and 
Co n (rf£ 5 d7 2 , S = 3/2)EB ra (dE 5 ,S = 1/2), respectively 
These are depicted in Fig.g as LS0 and HS0 states. 
Figure schematically represents the most plausible el- 



ementary processes in the cycle of the photo-induced 
structural change. The LH0 and HS0 states are con- 
verted to the intermediate states LSI and HS1, respec- 
tively, by photo-induced charge transfer (CT) between 
iron and cobalt atoms, and then to the final states HS0 
and LS0 by intersystem crossing due to spin-orbit cou- 
pling at cobalt sites. The major structural difference be- 
tween the LS and HS states lies in the Co-N bond length 
ri(Co-N), which is longer by about 0.2A in HS than in LS 
p3| . This implies that the transition between LS and HS 
states requires a volume change, which may be a major 
origin of a large energy barrier between them, as in the 
case of a spin-crossover complex |l4|,[l5| ■ 

Based on these observations, we adopt a cluster model 
consisting of a cobalt ion and surrounding cyano ligands 
with a variable bond length between Co and N. Previ- 
ously we carried out first-principles band structure calcu- 
lations 0. It was found that the relevant d-bands have 
quite narrow widths compared with the level spacings, 
implying that the d-electrons are fairly localized with 
weak charge transfer character in this material. This al- 
lows us to use the cluster approach in the present paper. 
We have checked that cluster calculations produce results 
that are consistent with band structure calculations . 

To take into account the partial substitution of lig- 
ands by water molecules, wc consider the following three 
clusters with various numbers Ny/ of H2O substitutions: 
(a) Co(NC) 6 (N w =0), (b) Co(NC) 5 -H 2 (iV w =l), and 
(c) Co(NC) 4 ■ 2H 2 (iV w =2). Bond lengths other than 
Co-N are fixed at observed values, e.g., d(C-N)=1.12A 
The potential energies of the four relevant states 
shown in Figj| are calculated for these clusters by using a 
quantum chemical ab-initio calculation program with the 
MIDI basis function set jL(|[l7). The multi-configuration 
self-consistent field (MCSCF) method is used with an 
active space limited to five orbitals mainly composed of 
Co-d atomic orbitals. The spin orbit coupling was not 
included to obtain the energy curve for each spin state. 

As for the treatment of surrounding atoms outside the 
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cluster, we use a lattice of atomic point charges evaluated 
from our previous band calculation H). The Fe-C bond 
length is fixed at the observed value 1.91 A |l3|. The 
environment lattice contains one or two iron vacancies 
situated just outside the water molecules of the clusters 
(b) and (c). The energies of the four states are calculated 
for each value of d(Co-N) common to all the Co-N bonds 
in the cluster as well as in the environmental lattice. This 
choice diminishes the interfacial strain energy due to lo- 
cal volume change, being a suitable way to extract the 
tendency of the local cluster. It should be noted that the 
real relaxation dynamics is determined not only by the 
cluster potentials but also by (mainly elastic) interactions 
between clusters. 

Since each Fe has full six coordination of CN with 
the fixed Fe-C distance, we assume that the energy dif- 
ference between the two states Fe 11 (cfe 6 ,5 = 0) and 
Fe m (cfe 5 ,5 = 1/2) is a constant A irrespective of Co 
states. Then the energy difference between the LSO 
and HSO states, for example, of a Co-Fe pair is writ- 
ten in an additive form: Ai?LSo-HSO=-£'[Co III (de 6 , 5 = 
0)] -E[Co u {de 5 dj 2 , S = 3/2)] - A. We have checked the 
appropriateness of this approximation by a supplemen- 
tary calculation for a two- metal cluster Co(NC)s-(NC)- 
Fe(CN)5. In the following we use A=0.4eV. This choice 
ensures that the ground state is in LS state at x=0.3 (j]] 
and in HS state at x=0.5 [Q as observed experimentally. 

Figure || shows calculated cluster potentials as func- 
tions of ci(Co-N) for Nw = and 1. We see that the en- 
ergies of the charge-transferred states (LSI and HSO) are 
lowered in the cluster of Nw = 1 compared with the case 
of iV\y = 0. They become further stabilized in the cluster 
of JVw = 2 (not shown here) . This implies that the lig- 
and substitution enhances the electron affinity of Co. On 
the other hand, the energy differences between different 
spin states (between LSO and HS1 and between HSO and 
LSI) are rather insensitive to Nyj, indicating that the 
ligand substitution does not affect the Hund coupling 
very much. From the calculated cluster potentials we 
obtain the Frank-Condon excitation energies of the CT 
excitations (LSO^LSl and HSO^HSl). The results are 
summarized in Table [| together with the observed main 
absorption peaks and the excitation energies responsible 
for the structural change. The calculated excitation en- 
ergy ~ 2.3 eV of LS0->LH1 in Co(iV w =l) reasonably 
corresponds to the observed absorption energy 2.4 eV in- 
ducing the magnetization. This implies that the photoin- 
duced LS— >HS transition is triggered by CT excitation 
mainly between a Co atom with N-w=l and a nearby Fe 
atom. Note that the largest fraction (37%) of Co is in the 
configuration of iVw=l in Ko.4Coi.3Fe(CN)g, if iron va- 
cancies are randomly distributed in the crystal. The Co 
atoms with Ny/=0 (21% fraction) are disadvantageous 
for this process because of the large energy difference be- 
tween the metastable HSO state and the stable LSO state 
with a small potential barrier for HSO, as one can see in 
Fig. |(a). 

In the case of the reverse (HS^LS) transition, the ex- 



perimental photon energy (0.9 eV) required for demag- 
netization is much lower than the energy of the observed 
broad absorption around 2.3 eV. The latter may corre- 
spond to the calculated CT excitations at ~1.4 eV in 
Co(iV\v=l) and at ~2.6 eV in Co(iVw=2), whereas the 
former is assigned to the excitation energy of Co(iVw=0) 
at ^0.8 eV. This implies that the HS— >LS transition 
is triggered mainly by the excitation of Co atoms with 
Nyf=0, in contrast to the forward transition. This is 
reasonable, because the final LSO state is much lower in 
energy than the initial HSO state in the case of Nw=0 
shown in Fig. ||(a). 

The suggested main relaxation paths in the reversible 
transitions are indicated by arrows in Fig. § The path 
for the LS— >HS transition is shown in the case of iV\y=l, 
and the path for the reverse transition is indicated in 
the case of iVw=0. The Frank-Condon photo-excitations 
induce CT between cobalt and iron cations. The inter- 
mediate excited states are subsequently transformed into 
local meta-stable states via intersystem crossing due to 
the spin-orbit interaction at cobalt site. The magnitude 
of the spin orbit coupling has been estimated as about 
10 meV |l9| . This is much smaller than that of the typ- 
ical energy scale of the potential diagram in Fig. |. On 
the other hand, it corresponds to a time scale less than 
lps, which is much faster than the time scale of the pho- 
toinduced phase transition. Experimentally, intersystem 
crossing has been observed to occur with large probabil- 



ities in similar cobalt complexes |20 



The most important conclusion from our results is that 
the forward and the backward changes are initiated from 
spatially different local regions in the crystal. The dif- 
ference lies in the cluster potential for Co electrons due 
to the difference in the number of ligand substitutions: 
A r w=0 or 1. The large concentration of Fe vacancies 
in this crystal, which might otherwise be thought to be 
a shortcoming of this material, is an essential prereq- 
uisite for the reversible transformations. The present 
model also implies that the photoinduced change must 
be a cooperative phenomenon [ p~5|j2l] — 12^] , because the 
structural change of the entire crystal is induced by ex- 
citations in a fraction of the crystal. To confirm our 
model, a systematic experimental study with varying x in 
Ki_2a;Coi+ a ;Fe(CN)6 is desired especially. For instance, 
the oscillator strength of the 3.9eV absorption peak in 
the LS state will be much more enhanced for smaller x 
because of an increased fraction of Co(-/V\y=0). 

Based on the present results, we are able to specu- 
late how the initial relaxation process leads to the phase 
transition in the macroscopic scale. A local relaxed CT 
state (exciton) after the primary process is still in an 
excited state with a finite life time in the crystal, be- 
cause it involves a strain energy due to lattice mismatch 
or frustration. The strain energy is released when the 
surrounding lattice undergoes a global change as an inte- 
gration of all the local changes. There should be a crit- 
ical concentration, Rq, of relaxed excitons, above which 
the fraction of the new phase increases drastically in the 
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crystal. Such behavior was observed in a similar photo- 
induced phase transition of a spin-crossover complex [fl5| . 
The exciton concentration R satisfies P PI (1 — R) = R/t 
in the stationary condition, where P PI and r represents 
the creation rate and the lifetime of exciton, respec- 
tively. The critical excitation rate Pq 1 is then expressed 
as P PI = Rc/[t(1— Rc)]- Therefore, the phase transition 
occurs more easily with a smaller excitation light inten- 
sity for a longer exciton lifetime. In the present system, 
the long life time of the relaxed exciton is guaranteed by 
its local stability as shown in Fig. [| 

The local mechanism discussed here indicates a general 
strategy to realize a reversible photoinduced phase tran- 
sition by the use of a mixed crystal of a and b molecules. 
(Here a 'molecule' means just a structural unit in a wide 
sense.) Suppose both the a and b molecules are bistable 
between the two types of state A and B, but they are dif- 
ferent in that the a molecule prefers the A-type state and 
the b molecule prefers the B-type state. (In the Prussian 
blue analogues the clusters with iVw=l and play the 
roles of the a and b molecules.) In the A phase of the crys- 
tal where all the molecules are in the A-type states, the b 
molecules are metastable locally. By selectively exciting b 
molecules by light, their states are readily changed to the 
locally stable B-type states, triggering the phase transi- 
tion of the entire crystal to the B phase where all the 
molecules are in the B-type states. In the B phase, the 
story is the same with the interchanged roles of the a and 
the b molecules: Photo-excitation of locally metastable a 
molecules triggers the phase transition to the A phase. 

In summary, our ab initio calculations have elucidated 
that the photoinduced forward and backward phase tran- 
sitions are initiated by excitation of different local regions 
by different photons in Prussian blue analogues. Espe- 
cially, the existence of a large concentration of vacancies, 
which might be thought as a shortcoming of this mate- 
rial, turned out to be essential for the bidirectional phase 
transition. We have also proposed that this mechanism 
can be generalized as an effective strategy to design a 
new material for bidirectional photoinduced switching. 
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TABLE I. Frank-Condon excitation energies of the LS and 
HS states evaluated from the cluster potentials (as shown in 
Fig. 3) for various Aw, compared with experiments of the 
absorption peak energies (underlined) and the excitation en- 
ergies of the phase transitions (in bold). 
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FIG. 1. Crystal structure of Ki_ 2;l! Coi +x Fe(CN) 6 : (a) 
x — 0; (b) x =^ 0. Black, white, and gray circles denote Fe, 
Co and CN, respectively. In (b), there is an vacancy at an 
iron site surrounded by six water molecules substituting cyano 
anions. 
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FIG. 2. Schematic local spin configurations of a Co-Fe 
pair. LSO, LSI, HSO, and HS1 represent the non-magnetic 
ground state, the paramagnetic excited state, the mag- 
netic meta-stable state, and the magnetic excited state, re- 
spectively. The magnetization LSO^HSO occurs through 
charge transfer excitation to LSI by photon hv and subse- 
quent intersystem crossing (I.C.), whereas the demagnetiza- 
tion HSO^LSO occurs similarly via HS1. 
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FIG. 3. Calculated potentials of the cobalt-centered clus- 
ters with (a) iV w =0 and (b) N w =l. The LSO, LSI, HSO, HS1 
states are shown by solid, broken-dotted, dotted, and broken 
lines, respectively. Possible primary relaxation paths for the 
LSO^HSO and HSO— >LS0 transitions are schematically indi- 
cated by arrows in (b) and (a), respectively (see text). 



•5 



